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ABSTRACT 

We report the discovery of a galaxy cluster serendipitously detected as an extended X-ray source 
in an offset observation of the group NGC 5044. The cluster redshift, z=0.281, determined from the 
optical spectrum of the brightest cluster galaxy, agrees with that inferred from the X-ray spectrum 
using the Fe Ka complex of the hot ICM (z = 0.27 ± 0.01). Based on the 50 ks XMM observation, 
we find that within a radius of 383 kpc the cluster has an unabsorbed X-ray flux, fx (0.5 — 2keV) = 
(3.34lg;5|) x 10~ 13 ergs cm" 2 s~\ a bolometric X-ray luminosity, L x = (2.21+§;?£) x 10 44 ergs s~\ 
kT = 3.57 ± 0.12keV, and metallicity, 0.60 ± 0.09 Zq. The cluster obeys the scaling relations for Lx 
and T observed at intermediate redshift. The mass derived from an isothermal NFW model fit is, 
M v i t = 3.89 ±0.35 x 10 14 M Q , with a concentration parameter, c = 6.7 ±0.4, consistent with the range 
of values expected in the concordance cosmological model for relaxed clusters. The optical properties 
suggest this could be a "fossil cluster" . 

Subject headings: galaxies: clusters: general — X-rays: general 
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1. INTRODUCTION 

Galaxy groups and low mass clusters (T < 4 keV) are 
starting to be detected and analyzed in detail in the X- 
ray band at intermediate redshift 0.2 < z < 0.6 (e.g., 
IWillis et all 120051 : Ueltema et aD 120061 ). These objects 
are more likely to display the effects of non-gravitational 
energy injection into the intracluster m edium (ICM ) than 
hotter, more massive, clusters (e.g., |Voit| [20051). The 
study of extended X-ray objects over a broad tempera- 
ture range at z > 0.2 will provide important insight into 
the evolution of their hot gas and the deviation of X-ray 
scaling relations from simple, self-similar expectations. 
Studies of objects in the redshift range 0.2 < z < 0.6 
with kT ^2 — 3 keV already suggest that they are less 
dynamically evolved t han their counterparts at z = 
(jMulchaev et alJl2006f >. 

The combination of surveys specifi cally designed to de - 
tect clusters (e.g. the XMM-LSS. iPierre et all |2004[ ). 
with serendipitous observations, which can make use of 
deeper exposures, will further advance the knowledge in 
this field. Here we present the discovery of a cluster at 
the redshift of z = 0.281, serendipitously observed dur- 
ing an offset observation of the nearby group NGC 5044. 
All distance-dependent quantities have been computed 
assuming Ho— 70 km s _1 Mpc -1 , fi m = 0.3 and Sl\= 
0.7. At the redshift of z = 0.281, 1' corresponds to 255 
kpc. All the errors quoted are at the 68% confidence 
limit. 

2. X-RAY ANALYSIS 

The object XMMU J131359.7-162735 was detected in 
an offset observation of the group NGC 5044 (obsID 
0301290101) performed on 6 January 2006. The source 
(see FigfT]) is located at 4.7' off-axis, and it is clearly ex- 
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tended on a ~ 2 arc-minute scale. The X-ray centroid 
of XMMU J131359. 7-162735 in equatorial coordinates is 
aj2ooo.o = 13 h 13 m 59.7 s , «5 J20 oo.o = -16°27 m 35 s . 

The data were reduced with SAS v7.0.0 using the tasks 
emchain and epchain. We considered only event patterns 
0-12 for MOS and for pn. We cleaned the data using 
the standard procedures: bright pixels and hot columns 
were removed by applying the expression FLAG == 0, 
and we corrected for pn out-of-time events. Periods of 
high background due to soft protons were filtered as in 
iGastaldello et al. (2006); the observation was very mildly 
affected by flares leading to only ~ 5 ks of lost data, with 
net exposures of 52, 50 and 45 ks respectively for MOS1, 
MOS2 and pn. 

For each detector we created images in the 0.5-2 keV 
band with point sources masked out using circular re- 
gions of 25" radius centered at each source position. 
The images have been exposure-corrected, and a radial 
surface brightness profile was extracted from a circu- 
lar region of 6' radius positioned at the cluster cen- 
troid. We account for the X-ray background and the 
emission of the NGC 5044 group by including a con- 
stant background component. The data were grouped 
to have at least 20 counts per bin in order to apply 
the x 2 statistic. The fitted model is folded with the 
XMMPSF at energy of 1 keV. The joint best-fit model 
(|Cavaliere fe Fusco-Femianolll978f) has a core radius of 
r c = 66 ± 6 kpc (15.5"±1.4") and = 0.54 ± 0.01 for 
X 2 /d.o.f. = 465/273. Although formally unacceptable, 
the surface brightness model well-describes the overall 
shape of the profile, as evinced by the low fractional 
residuals (see Fig(2]) . Fits to the profiles of the indi- 
vidual detectors give consistent results within la of the 
combined-fit result. If we use larger extraction regions, 
the emission by NGC 5044 can no longer be treated as 
constant: if we account for the emission of NGC 5044 by 
adding a second model at an offset of 20.7' (the distance 
between the source and the e mission peak of NG C 5044) 
with r c = 48" and = 0.53 (jBuote et al.ll2003l ). the re- 
sulting parameters for the source are consistent within 
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Fig. 1.— Left: Exposure corrected 0.5-2.0 keV combined MOS1 and MOS 2 X-ray image of the observed offset field of NGC 5044. The 
cluster XMMU J131359. 7— 162735 is clearly visible as the extended source. Right: The V band image with the smoothed X-ray contours 
overlaid and logarithmically spaced between 2 and 37 cts/pixel. Coordinates are J2000. The red circle shows the FWHM (6.7") of the 
XMM MOS PSF at the off-axis angle of the source. 



the la errors. 

We extracted spectra for each detector from a 1.5' 
circular region positioned at the centroid of the emis- 
sion, chosen to maximize the S/N over the background. 
(Choosing an extraction radius of 2' does not change 
the spectral parameters quoted below). Redistribution 
matrix files (RMFs) and ancillary response files (ARFs) 
were generated using the SAS tasks rmfgen and arfgen, 
the latter in extended source mode. Appropriate flux- 
weighting was performed for RMFs, using our own ded- 
icated software, and for ARFs, using exposure-corrected 
images of the source as detector maps (with pixel size 
of 1', the minimum scale modeled by arfgen) to sam- 
ple the_j^anatioji_JnenTission, following the prescription 
of ISaxton fc Siddiquil (|2002). The background was es- 
timated locally using spectra extracted from a 2'-3' an- 
nular region positioned at the centroid of the emission; 
using a 3'-4' annulus gives results consistent within their 
la errors. We checked that vignetting corrections for 
the sky background components do not affect the results 
by exploring the complementary approach of a complete 
modeling of the background compone nts (including the 
therm al emission of NGC 5044) as in Gastaldell o et all 
(2006): the spectral results obtained are consistent with 
the ones quoted below within their la errors. 

The spectra from the three detectors were re-binned 
to ensure a signal-to-noise ratio of at least 3 and a min- 
imum 20 counts per bin, and they were jointly fitted 
with an APEC thermal plasma mod ified by Galactic ab- 
sorption (Dick ev fc Lockmanl |1990T ) . T he spectral fit- 
ting was performed with Xspec (verll.3.1. lArnaudlll996t ) 
a nd quoted metallicit i es are relative to the abundances 
of iGrevesse fc Sauv al (1998). We display the spectra in 
FigEl The source photons correspond to about 77% 
of the total events (~ 4400 counts in each MOS and 
~ 7700 in the pn). An emission line is clearly visible 
corresponding to the Fe Ken complex emitted at a red- 
shift of z — 0.27±0.01. The redshift measured is in close 
agreement with the redshift of z — 0.2814 determined op- 
tically from the spectrum of the brightest cluster galaxy 
(see §3j) - Fixing the redshift value to the optical deter- 
mination gives as best fit parameters, T — 3.57 ± 0.12 
keV, and, Z = 0.60 ± 0.09 Z @ , for x 2 / d -°- f - = 486/463. 

Using the best-fit model, the unabsorbed flux within 
the aperture of radius 1.5' (383 kpc) is 3.34±°;^ x 10~ 13 
ergs cm -2 s~ 4 in the 0.5-2 keV band. This corresponds 



to an unabsorbed luminosity of 7.571q' 4 q X 10 43 ergs s 



in the 0.5-2 keV band and to a bolometric (0.01-100 keV) 



luminosity of 2. 2llg' 3 g x 10 44 ergs s _1 . The quoted errors 
on flux and luminosity are obtained by Xspec using a 
Montecarlo procedure. 

To investigate possible spatial variation in the spectral 
parameters of the cluster, we extracted two annular re- 
gions of radii 0'-0.5' and 0.5'-1.5'. The derived spectral 
parameters are: T = 3.64 ± 0.14 keV and Z = 0.66 ±0.13 
Zq with x 2 /d.o.f. = 280/231 for the inner annulus; 
T = 3.48±£^ keV and Z = 0.52 ±0.13 Z with x 2 /d.o.f. 
= 280/271 for the outer annulus. The width of the bins 
has been chosen in order to avoid bias in the temperature 
measurements caused by scattered flux by the PSF ( 80% 
encircled energy fraction radius is 25" at 1.5 keV and at 
the off- axis angle of the source). The cluster is therefore 
consistent with being isothermal over the explored radial 
range. 

XMMU J131359. 7-162735 was also serendipitously de- 
tected in the ROSATPSPC pointed observation of NGC 
5044 (1WGA J1313. 9-1627: White et al1fl99l . though 
it was very close to the detector gaps. This is likely 
the reason why with RO SAT the source appeared not to 
be extended and had a lower flux (1.75 ± 0.14 x 10~ 13 
ergs cm -2 s _1 in the 0.2-2 keV band). 

The cluster has regular X-ray isophotes and is cen- 
tered on a dominant early type galaxy (see SJ2]). These 
characteristics suggest the cluster is relaxed and that hy- 
drostatic equilibrium is a good approximation. We cal- 
culated the total mass profile using two different mod- 
els. First, we used the best-fit model (see 32) f° r 
which the gas density and total mass profi les can be ex - 
pressed by simple analytical formula (e.g.. lEttorill2000j ). 
We evaluated r^o &s the radius at which the density is 
500 times the the critical density and the virial radius 
as the radius a t which the density corr esponds to A v j r , 
as obtained bv iBrvan fc Normanl (|1998h for the concor- 
dance cosmological model used in this paper. To evalu- 
ate the errors on the estimated quantities, we repeat the 
measurements after 10000 random selections of a tem- 
perature and parameters of the surface brightness pro- 
file, which were drawn from Gaussian distributions with 
mean and variance in accordance with the best- fit results. 
For A = 500 we obtained, M 500 = (1.60 ± 0.09) x 10 14 
Mq within rsoo = 750 ± 14 kpc; the virial mass is, 



XMM discovery of a cluster of galaxies at z = 0.28 



3 




Radius (arcmin) 

Fig. 2. — Surface brightness profile of the X-ray emission of 
XMMUJ131359. 7-162735. Data from MOS1, MOS2 and pn arc 
plotted in green, red and black respectively. The best fit beta model 
and ratio of data over the model are also shown. The background 
level for PN and for one of the MOS are shown as black and green 
dashed lines respectively. 

M vir = (3.24 ± 0.18) x 10 14 M Q , within the virial ra- 
dius TVir = 1511 ± 29 kpc. Secondly, we fit the sur- 
face brightness p rofile with an isothermal NFW model 
(jSuto et all IT998). We obtain a concentration parame- 
ter, c = 6.7 ± 0.4, virial radius r vir = 1607^42 kpc, and 
virial mass M vir = (3.89±0.35) x 10 14 M Q , with x 2 /d.o.f. 
= 424/273. The mass determinations using the two dif- 
ferent models are in agreement within the ler errors. 

3. OPTICAL FOLLOW-UP OBSERVATIONS 

A I^-band and /-band image of the field were taken 
with the direct imaging CCD on the 2.5-m du Pont 
Telescope at Las Campanas Observatory on the night 
of 27 January 2006. In both wavebands, 2x300 s expo- 
sures were taken in clear, photometric conditions. Pho- 
tometric zero-points and extinction corrections were de- 
ter mined by observ ing three standard stars chosen from 
the lLandoltl ([1992) catalog. The images have been pro- 
cessed and cleaned of CCD defects and cosmic rays. The 
y-band image is shown with X-ray contours superposed 
in FigQ] The image reveals an over-density of galaxies 
and, in particular, a bright dominant galaxy at the same 
location as the X-ray centroid. 

Photometry was pe rformed using SExtractor 
(jBertin fe Arnoutsl Il996f) with automatic aperture 
magnitudes. We plot in the top panel of Fig[4] the 
color-magnitude diagram of the galaxies detected within 
3.1' radius of the central galaxy (corresponding to 803 
kpc, half of the estimated virial radius): the horizontal 
dashed line corresponds to the expected V — I color of a 
passively evolving elliptical galaxy at z ~ 0.28, which is 
V — I ~ 1.6 for a formation redshift, z > 5 ([Nelson et al.l 
2002). The difference in magnitude between the central 
galaxy (mi — 16.4) and the likely cluster members 
is very close t o the 2 magnitud es of the definition of 
fossil systems ([Jones et al.1 120031 ). We have excluded 
three objects (triangles in FigHJ) because they are likely 
members of the NGC 5044 group. A further two objects, 
which are within 2 magnitudes of the BCG, are very far 
from the cluster center (2.6', 664 kpc, and 2.5', 638 kpc, 
crosses in FigfJJ. A possible member is 1.9 magnitudes 
fainter than the BCG at a distance of 1.8' (diamond in 
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Fig. 3. — X-ray spectrum of the source taken from a 1.5' aperture 
positioned at the centroid of the emission. Data from MOS1, MOS2 
and pn are plotted in green, red and black respectively. The best 
fit model and residuals are also shown. 

FigiD. 

We obtained an optical spectrum of the central galaxy 
on the night of 7 March 20 06 UT, using the ESI spectro- 
graph ([Sheinis et al.l 1200^ 1 at the Keck-II telescope. A 
0'.'75 slit width was used, giving a spectral resolution of 
~ 50 km s _1 (FWHM) and a pixel scale of 11.4 km s _1 
pixel -1 . The galaxy was observed at airmass 1.25 for 900 
s in partly cloudy conditions, and the slit was oriented at 
the parallactic angle. The spectrum was extracted with a 
1" extraction width, and wavelength calibrated and flux- 
calibrated using an observation of the standard star Feige 
34. The individual echelle orders were combined into a 
single spectrum. The galaxy redshift was measured by 
performing a direct fit of a K3III stellar spectrum (ob- 
served on the same night with ESI), velocity-broadened 
by convolution with a Gaussian kernel. The best-fitting 
model yields a recession velocity of 84360 ± 19 km s -1 , 
or z = 0.2814 ± 0.0001, in close agreement with the red- 
shift estimated from the X-ray observations. The ESI 
spectrum is shown in the bottom panel of FigfJ] 

4. DISCUSSION 

This serendipitously discovered cluster obeys the scal- 
ing relations for Lx and T derived at intermediate red- 
shift. The aperture of 383 kpc used for spectroscopy en- 
closes 78.3% of the flux within rsoo, assuming the cluster 
emission profile follows the [3 model derived within that 
aperture. The derived bolometric luminosity within f5oo 
is Lsoo = (2.82t°;g) x 10 44 ergs s" 1 . The derived L 500 is 
in good agreement with the L — T relation f ound for the 
sampl e of low redshift clusters (z < 0.09) by Mar kevitchl 
(|1998f) . the low-redshift groups (z < 0.03) in the GEMS 
sample ([Osmond fc Ponmanll2004T ). the six groups/poor 
clusters at intermediat e redshift (0-29 < z < 0.44) in 
the XMM LSS survey (| Willis et al.ll2005l ). and the six 
groups/poor clusters at i ntermediate redshift (0.23 < 
z < 0.59) in the sample of iJeltema et al.l (|2006f) . 

The optical appearance of the cluster suggests tanta- 
lizing evidence that this object may be one of the rare 
fossil clusters (3%-6% of the total number according to 
Milosavlic vic et al"1l2006l ). More detailed optical imaging 
and spectroscopy are required to investigate this possibil- 
ity. The derived total magnitude for the BCG, corrected 
for Galactic extinctio n, as well as for the fc-correction 
(|Fukugita et al.l I1995D . is Mi = -24.78, which corre- 
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Fig. 4. — Top: Color-magnitude diagram of the galaxies within 
2.8' of XMMU J131359. 7-162735. The horizontal dashed line in- 
dicates the predicted color of a z = 0.28 cluster elliptical. The 
star indicates the BCG. Triangles indicate likely dwarf members 
of NGC 5044 (one of the m is the object #26 in the catalog of 
Ferguson & Sandagc 1990). Crosses indicate objects very far from 
the BCG (see text). The diamond indicates the only other ob- 
ject within 2 magnitudes of the BGC. Bottom: A portion of the 
Keck ESI spectrum, binned to 2 A/pixel. Over-plotted in red is 
the elliptical galaxy template spectrum from the spectral atlas of 
IKinnev et all JT996T) . 



sponds to Mb = —22.38 and Mr = — 24.05, adopting 
color s B-Ij = 2 AO and Rj - Ij = 0.73 (jFukugita et all 
QUI). The BCG of XMMU J131359. 7-162735 is among 
the mo st luminous in the s ample of BCGs in fossil sys- 
tems of iKhosroshahi et alJ (|2006l ). The measured values 
of Lb of the BCG and Lx pu t the object very clos e to 
the locus of baryonic closure (|Mathews et alJ 120051 ) . as 
expected for fossil systems. 

The measure d c, multiplied by t he expected depen- 
dence of 1 + z (|Bullock et al .11200 If ), is consistent with 
a ACDM model with a$ = 0.9 indicative of a relaxed, 
early forming system in agreeme nt with the theoret- 
ical prediction for fo ssil systems (jZentner et alJ 120051 : 
iD'Onghia et al.ll2005l ) and with the observat ional results 
of relaxed, low- 2 clusters (jBuote et al.ll2006l) . 

5. CONCLUSIONS 

The discovery of XMMU J131359.7-162735 reported in 
this letter represents the deepest observation among the 
intermediate redsh i ft systems discussed in the literature 
(jWillis et all 120051 : I Jeltema et al] 120061) . and it already 
provides an interesting description of its fundamental 
parameters and dynamical state. Future optical spec- 
troscopy and dedicated X-ray follow-up will allow de- 
tailed spectral properties to be measured. 
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